Instantaneous Measurement of field quadrature moments and entanglement 
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We present a method of measuring expectation values of quadrature moments of a multimode field 
through two- level probe "homodyning" . Our approach is based on an integral transform formalism 
of measurable probe observables, where analytically derived kernels unravel efficiently the required 
field information at zero interaction time, minimizing decoherence effects. The proposed scheme 
is suitable for fields that, while inaccessible to a direct measurement, enjoy one and two-photon 
Jaynes-Cummings interactions with a two-level probe, like spin, phonon, or cavity fields. Available 
data from previous experiments are used to confirm our predictions. 
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State reconstruction of a bosonic field is an impor- 
tant issue in fundamentals of quantum physics that has 
been studied extensively, both theoretically and experi- 
mentally, in the last two decades Q- Its main concern 
is to measure either the density matrix of an unknown 
quantum field state or, equivalently, any of its phase- 
space representations. Among them, Wigner function 
reconstructions seem to be the most promising avenue 
for measuring completely, for example, an intracavity mi- 
crowave field or the motion of a trapped ion Q . In the 
case of a propagating field, usually accessible to direct 
measurement, homodyne techniques are currently used 
in the lab |4j. Typically, a complete state reconstruction 
with standard techniques demand great experimental ef- 
forts, is strongly affected by decoherence mechanisms, 
and, frequently, the obtained information exceeds our 
requirements. In those cases, techniques for extracting 
efficiently the required partial information are most wel- 
comed and even necessary. The problem is even harder 
when the field is not directly accessible and a quantum 
probe has to be used for the purposes of an indirect mea- 
surement 0. Therefore, the following question arises: 
how to derive accurately the expectation value of a field 
observable, through an efficient measure of a probe, with 
minimal resources and with an outcome that is minimally 
affected by decoherence mechanisms? In this article, we 
answer this question for the case of a multimode bosonic 
field, interacting with a two-level probe, by means of a 
practical integral transform method. These conditions 
are naturally fulfilled by several physical systems, like a 
cavity field interacting with two-level atoms, the motion 
of an ion interacting, through laser coupling, with two of 
its internal levels, or even several spins, in a mode ap- 
proximation, interacting with a single spin, like in NMR 
or quantum dot systems. 

We consider a general picture in which an inaccessible 
bosonic field is measured through an interacting probe, 
following the interaction Hamiltonian 



(1) 



where pi and fj are probe and field operators, respec- 
tively, and gi j are coupling strengths. We postulate the 



existence of an analytical kernel k(t) such that 

/•OO 

(F)= K(T)(P)(r)dT, (2) 



where F and P are operators associated with field and 
probe observables, respectively, and r is the dimension- 
less probe-field interaction time. Later, it will be clear 
why it is possible to include a negative axis for r in the 
integration limits of Eq. |J2J. (P)(t) will be replaced 
by experimental measured data and, for the method to 
be useful, we should be able to formally invert the in- 
tegral transform and derive an analytical expression for 
k(t). We will show below that this inversion is possi- 
ble for several important field observables, like quadra- 
ture moments of a multimode field, unravelling impor- 
tant information on squeezing and entanglement. Note 
that we aim at measuring efficiently partial field informa- 
tion without the requirement of full state reconstruction, 
even though integral techniques can also provide us with 
complete Wigner function reconstructions @- 

We consider a two-level probe interacting with a single- 
mode field through a resonant Jaynes-Cummings (JC) 
Hamiltonian, in the interaction picture, 



H 



JC 



hg(a' i a + era'), 



(3) 



where g is a coupling strength, {<r, er^} are lowering and 
raising probe operators, and {a, a 1 } are annihilation and 
creation field operators. We assume that, given that the 
initial probe-field density operator is p m (0), we can mea- 
sure, after a dimensionless interaction time r = gt, the 
population of the excited probe level |e) 

J*» = Tr[p(r)\e){e\] = Tr[U(r)p in (0)U^ (r)\e){e\], (4) 

where \e)(e\ = a^a and U(t) = exp(— irHjc/hg) is the 
evolution operator. In Eq. and throughout this work, 
upper and lower indices stem from initial and measured 
probe states, respectively. We consider the initial state 
P + " = \+<p)i+<p\®Pf, where |±«A = (\g) ± e*\e))/y/2 are 



the eigenvectors of erf 



ith 



When 4> — 0, <j% turns into a Xl the usual spin- 1/2 
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Pauli operator. Replacing p + * in Eq. J2J gives 

. oo 

P^( T ) = sin(2rV n + l)(e'^p re ,tt+i - e~ l< V„ + i, n ) 
n=0 

_^ oo 

+ - ^(cos 2 (rVn~TT)p n ,„ + sin 2 (rVn + I)p n+ i, n+ i), 

(5) 

where the initial field state = ^ n m P«,m| n )( m l nas 
been written in terms of its matrix elements. Through 
the knowledge of P e *(t), we aim at measuring the field 
quadratures X^ = ^{ae~ l< ^ + a^e" 1 ^) and = X,p +7r /2 — 



i(f> 



ve i( P) with expectation values 
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= ? J2 V^TTfe* + e^V+i,™) (6) 

(7) 



n=0 

oo 



( Y 4>) = 7} Vrt + l(e i0 /9„ in+ i - e 

n=0 



We choose the kernel in Eq. as an odd function, 
k(—t) = —k(t), so that the integral of the second sum 
in Eq. (JSJ) vanishes, while the integral of the first sum 
should reproduce Eq. Q. In consequence, replacing op- 
erators P — > |e)(e| and F — > Y^ in Eq. J3J), the condition 
for this ansatz to be true is 



(8) 



where sin(2r-\/n + 1) has been rewritten in complex form 
with w n = 1\Jn + 1. The inverse Fourier transform of 
Eq. JHJ provides us with the kernel 



«(r) 



i 

2^ 



e-^ T u> n du> n = -S'(t), (9) 



where S'(t) is the first derivative of a delta function. Note 
that even if, physically, w n is a function of discrete ra's, 
it can be treated formally as continuous for the sake of 
the inverse transform. Then, Eq. J3J) can be written as 



{Y l p) = - S'(r)P+*(T)dr, 



yielding 



(10) 



(11) 



where the continuity of the first derivative of P e * (r) at 
t = has been considered. Similarly, 



<^> = <Y^_ 5 ) = -P e *- f (T) 



(12) 



T = 



Eqs. {TT} and JEJ) show that (X$) and (Y^) are fully 
contained in the first derivative, at r = 0, of the measured 



probe population, offering a remarkably simple way of ob- 
taining quadrature information. Note that knowing the 
first derivative at r = requires knowing the function in 
a vicinity. However, no necessity of full state reconstruc- 
tion or lengthy time integration over Rabi oscillations 
are needed, in contrast to standard methods. Needless 
to say, the influence of decohcrcncc is minimized. 

Induced by the structure of Eq. (J5J, and aiming at 
cancelling the population while keeping the off-diagonal 
(phase) information, we could find a similar result by 
subtracting rotated populations 



Pe*(T)-Pe*(T) 



- ^ sin(2r%/nTT)(e 40 p„ in+ i 



— id> 
e Pn+l, 



(13) 



Following a similar procedure as before, we can write 



(14) 



r=0 



This result has evident resemblance to the known tech- 
nique of field homodyning There, an unknown field 
is mixed in a 50-50 beam splitter with a local oscilla- 
tor, and the difference of field intensities (rate of photon 
clicks) at the output gives us quadrature information. 
Based on this similarity, the proposed method could be 
called after two- level probe " homodyning" . 

Xrf, and Yrf, happen to be relevant observables in a wide 
range of physical systems where current experiments en- 
joy probe rotations and JC-like interactions, like cavity 
QED (CQED), trapped ions, BEC, and different solid- 
state systems. In CQED, the quadrature information 
can be obtained by sending an excited atom through a 
Ramsey zone before crossing the cavity mode 7], and 
finally measuring the population of the excited state at 
the cavity output. For trapped ions, (X) and (Y) repre- 
sent, literally, expectation values of position and momen- 
tum operators, that will be obtained by measuring the 
internal state statistics, where the efficiency can reach 
~ 100%, after a JC-like sideband excitation In the 
case of solid-state devices, there are several systems en- 
joying two-level probes interacting through JC interac- 
tions with cavity, phonon or spin fields. It is noteworthy 
to mention that in all these examples a probe is needed 
due to the lack of a direct measurement. 

In the rest of this article, for the sake of simplicity, we 
will use the language of cavity QED, where a two-level 
atom probes an intracavity electromagnetic field. 

Another important field observable that can be ob- 
tained straightforwardly with a JC interaction is the 
mean photon number (n) = (a^a). Considering the ini- 
tial state p c = |e)(e| ® pf, we can derive the kernel 



S(r) = -8"(r) 



for measuring 



1 d 2 P c (r) 
K 1 2 



d 2 i 



(15) 



(16) 
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Note that measuring (n) does not require Ramsey zones 
for rotating the atom, as was the case before. Given the 
available experimental data, expression in Eq. H16[) is the 
only one that could be presently tested. For example, 
using the experimental data associated with the experi- 
ments at ENS, see Figs. 2 (A), (B) in Ref. @], we predict 
(n) w 0.14 and 0.81, respectively. These values are quite 
close to the ones obtained via integration or fitting long 
Rabi oscillations, 0.06 and 0.85, respectively. We made 
similar estimations for the experiments at NIST, obtain- 
ing (n) ss 1.6 and 3.1 for the experiments associated with 
Figs. 2 and 3 in Ref. 0, to be compared with 1.5 and 2.9, 
respectively. Clearly, our predictions could only be bet- 
ter if specific experiments are performed, aiming at first 
and second derivatives at very short interaction times. 

It is also possible to use these integral methods to mea- 
sure second-order quadrature moments, providing infor- 
mation about field quadrature squeezing and entangle- 
ment of a multimodc field. We will use a resonant two- 
photon JC Hamiltonian that reads 

H 2 jc = hg(a t a 2 + aa^ 2 ) (17) 

in the interaction picture. This nonlinear interaction has 
been realized experimentally in the context of microwave 
CQED [Tjj and trapped ions [j|. Our aim, here, is to 
measure expectation values of squared quadratures, 

TO - 4 + — 
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Pn,n+2 + e 2l ^Pn+2,n) 



(18) 



4 ' ~2 

OO 



~]TV(n + l)(n + 2)( e 2< < 



i -2id> 
Pn,n+2 + e Pn+2,: 



(19) 

with the help of Eq. (|17fl and the proposed integral trans- 
form techniques. Then, in a close analogy to Eq. (|13fl . 
now for a two-photon JC interaction, we can calculate 

OO 

P+*{t)-P-*{t) = % l Y J sin(2rV(n + l)(n + 2)) 

71 = 

x(e 2l V„ +2 -t-e- 24 V +2 ,„). (20) 

With the help of Eqs. (02), |Q3J, and JEJl, and by deriv- 
ing and using the corresponding kernels, we arrive at 
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1 ( dP^(r) 


dp-*{T 


dr 
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ld 2 P|(r) 
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T = 



r=0 



(21) 



(22) 



The quadrature variances (AX) 2 — (X 2 ) — (X) 2 and 
(AY) 2 — (Y 2 ) — (Y) 2 contain information about field 
squeezing and can be calculated straightforwardly by us- 
ing Eqs. (EJ), I2TJ, and 

It is noteworthy to say that it is not necessary to use 
a two-photon JC interaction for measuring second-order 
quadrature moments. For example, it is enough to use a 
two-atom probe interacting with the tested field through 
a single-photon JC, whose interaction Hamiltonian reads 



H l = hg[(a\ + a\)a + (a x + a 2 )a 



(23) 



where the subindexes are labelling probe atoms "1" and 
"2". We consider the Bell states \<f>g) = [\g1g2) + 
e ie \eie 2 )]/V2 and \(j> e ) = [\g x g 2 ) - e w \e x e 2 )\l^2 as two 
probe initial states, and in both cases we measure = 
[\91e2) + |ei3 2 )]/v / 2, obtaining 



2^ 

n=0 



V( n + l)(n + 2) 
2n + 3 



x sin' 



1 (V2rv / 2^T3)(e i Vr l +2 + e' l6 p n+2 , n ). (24) 



From this expression, and following similar steps to pre- 
vious derivations, it is possible to deduce 



TO 



(Y e 2 ) 



Ud 2 P;i(r) 

d 2 T 



d 2 Pll{r) 



lfP g (r) 



T = 
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U d 2 P*i(T) d 2 P^(r) 
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d 2 i 



T=0 



T = 
1 

4' 

T = 
1 

4 ' 



(25) 



(26) 



Note that the required Bell states and the measurement 
procedure have already been implemented in the lab in 
the case of CQED [13 and trapped ion [H setups. 

The formalism for measuring squeezing can be gener- 
alized to a two-mode field (or more), providing us with 
entanglement information. Accordingly, we define the 
two-mode quadratures as 



Xa = Xa 



1 2 



i=i 



(28) 



where j labels modes "1" and "2". The quantities (Xa) 
and (Ya,) can be easily calculated and, here, we will con- 
centrate on the second-order quadrature moments 

TO = TO) + TO)+2(A 01 A 02 ) (29) 
(Y 2 ) = (Y 2 1 ) + (Y 2 2 )+2(Y 4 , 1 Y (f , 2 ), (30) 
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In these expressions, we define (X t p 1 X ( f >2 ) = \(A) + \(B), 
(Y*M = ±(A)-±(B),with 

A = a^e^ 1 -^ +a l a\e^ 1 -^\ (31) 
B = alale-^ 1 ^ + aia 2 e lWl+ * 2) . (32) 

Single- mode quantities (Xfa ) , (Y ( j >i ) , ) and (Y£. ) , can 
be determined using two-level probes as it was shown 
above. Therefore, the main issue is to calculate the ex- 
pectation values of A and B, which describe correlations 
between modes 1 and 2. It has been shown, theoreti- 
cally and experimentally [l^, that the two-photon 
probe-field interaction Hamiltonian 



(33) 



can be engineered and controlled. If the probe is prepared 
initially in the superposition states, |+^) or |— with 
4> = 4>i — cf>2, we can calculate 



5 S 1 

n l ,1^2—0 



in ( 2gr \/n 2 (ni + 1)) 



J4> , 



Pni,rt2;ni + l,n2 — l ~t~ ^ ^ Pni + l,rt2 — l;ni ,712 ) 5 (34) 

from which we can derive 



(A) = - 
9 



'<aW_^ 

Similarly, and by using the Hamiltonian 
Hb = hg(a^ a[a 2 + (712102), 



T=0 



we can deduce 

(B) = 



i ( dP^jr) P- B (r) 



(35) 



(36) 



(37) 



In consequence, we can also estimate two-mode field vari- 
ances (AX) 2 and (AY) 2 in terms of measurable probe 
observables. Furthermore, by using the same approach, 
we can compute the variances of EPR-like operators 



y Cl 1 v 

aoAi — i r A2, 

\ci \ a 

v 02 1 V 
aoYi - -, — r — r 2 , 

c 2 a 



(38) 
(39) 



where ao, Ci and c 2 are constants. For example, it was 
shown in Ref. 14] that a two-mode Gaussian state p is 
separable if, and only if, ((Au) 2 ) p + ((Av) 2 ) p > Oq + I/oq. 

In summary, we have shown how expectation values of 
quadrature held operators can be measured by means of a 
two-level probe, helped by a practical integral transform 
method and without the necessity of full state reconstruc- 
tion. Surprisingly, all relevant information is contained in 
first and second derivatives of measurable probe observ- 
ables at interaction time r = 0, making unnecessary long 
range probe measurements and minimizing decoherence 
effects. Also, we showed that a similar technique allows 
to measure second-order quadrature moments and vari- 
ances, that is, squeezing and entanglement. These results 
allow us to conjecture the possibility of realizing full state 
reconstructing with "instantaneous" measurements that 
are robust to decoherence. 
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